Small adiabatic polaron with a long-range electron-phonon interaction 
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Two-site single electron system interacting with many vibrating ions of a lattice via a long-range 
(Frohlich) electron-phonon interaction is studied in the adiabatic regime. The renormalised hopping 
integral of small adiabatic Frohlich polarons is calculated and compared with the hopping integral 
of small adiabatic Holstein polarons. 
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If phonon frequencies are very low, the local lattice 
deformation due to the electron-phonon interaction can 
trap the electron. This self-trapping phenomenon was 
predicted by Landau Q]. It has been studied in greater 
detail by Pekar 0, Frohlich Feynman Devreese 
5] and other authors in the effective mass approximation, 
which leads to a so-called large polaron. 

When the electron-phonon coupling is relatively 
strong, A = E p /D > 1, all electrons in the Bloch band 
are "dressed" by phonons. In this regime the electron ki- 
netic energy, which is less than the half-bandwidth (D), 
is small compared with the potential energy due to a local 
lattice deformation, E p , caused by electrons themselves. 
Here the finite bandwidth is essential, and the effective 
mass approximation cannot be applied. The electron is 
called a small polaron in this regime. The main features 
of small polarons were understood by Tjablikov @, Ya- 
mashita and Kurosava , Sewell j|, Holstein HI and his 
school 0, [n|, Lang and Firsov |12^ . Eagles |13| . and 
others and described in several review papers and text- 
books [1 El H E EHIl. In particular, in his pio- 
neering study of the small polaron dynamics Holstein 9] 
introduced a two-site molecular model of a single electron 
coupled with local (molecular vibrations). He derived a 
renormalised (polaron) hopping integral for two extreme 
cases, the nonadibatic, lu > t and adiabatic, uj <C t, where 
t is the bare (unrenormalised) hopping integral, and to is 
the characteristic phonon energy. An exponential reduc- 
tion of the bandwidth at large values of A is a hallmark 
of small polarons. The self-trapping is never "complete" 
even in the strong-coupling regime, that is any polaron 
can tunnel through the lattice. Only in the extreme adi- 
abatic limit, when the phonon frequencies tend to zero, 
the self-trapping is complete, and the polaron motion is 
no longer translationally continuous. 

In the last years quite a few numerical and advanced 
analytical studies confirmed these and subsequent 03 
analytical results ( see, for exa mple, lip. Efl Ell E^ . Ea 

HlMElElEllMlMIIIllil iSOllfr At the 

same time polarons were experimentally recognised as 
quasiparticles in novel materials, in particular, in su- 
perconducting cuprates and colossal magnetoresistance 
manganites [36| . 

Small polarons are very heavy due a local character 
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FIG. 1: Electron interacting with infinite chain of ions via 
long-range interaction performing a hop between two sites 1 
and 2. 



of the electron-phonon interaction in the Holstein model, 
so that their huge effective mass created some prejudices 
with respect to any relevance of polarons to real oxides. 
However, it was pointed out that small polarons could 
be quite mobile if one takes into account a more realistic 
long-range interaction with phonons in ionic solids |3j]. 
Indeed the Monte-Carlo [33 and other calculations [33, 
|39| proved that small polarons with a long-range Frohlich 
interaction are a few orders of magnitude lighter than 
the small Holstein polarons (SHP) with the same binding 
energy. All numerical and analytical results for these, so- 
called small Frohlich polarons (SFP) 33], were obtained 
in the nonadiabatic or near-nonadibatic regimes. Here 
we extend these studies to the adiabatic region. 

We consider an extended Holstein model of an elec- 
tron hopping between two sites, but interacting with 
all surrounding ions of the lattice via a long-range 
electron-phonon coupling, like in a simple case of a one- 
dimensional ionic chain, vibrating in the direction, per- 
pendicular to the chain, as shown in Fig.l. 

For simplicity, sites 1 and 2 are not vibrating, but the 
interaction with their vibrations can be easily included in 
our model. The model mimics high — T c cuprates, where 
in-plane (C1LO2) carriers are strongly coupled with the c- 
axis polarized vibrations of apical oxygen ions |40| . We 
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derive an analytical expression for the renormalised hop- 
ping integral and compare it with SHP in the adiabatic 
limit. 

The Hamiltonian of the model is 



H = H, 



ph 



ph 



where 



H ph - [~ 



d 2 



MuPu 2 
2Mdu 2 m + 2 

is the Hamiltonian of vibrating ions, 

H e -ph = Y y^qci/ m (i)t< 

i=l,2 m 



(1) 



(2) 



(3) 



describes interaction between the electron and the ions, 
and H e ~ — t(c\c 2 + H.c.) is the electron hopping energy. 
Here u m is the displacement and / m (i) is an interacting 
force between electron on site i and the m-th ion. M is 
the mass of vibrating ions and lo is their frequency. 

The wave function is a superposition of two normalized 
Wannier functions H^r), localized on the left and the 
right sites, 



* = ai (u m )W(r - 1) + a 2 (u m )W{r - 2) 



(4) 



The Schrodinger equation is reduced to two coupled 
second order differential equations with respect to the 
infinite number of vibrational coordinates u m 



E - J2( H ph ~ /m(l)Wm) 



ai(u m ) = ta 2 (um) (5) 



a2(«m) = tax(u m ) (6) 



Let us first discuss the nonadiabtic regime. Following 
Holstein [t| we apply a perturbation approach with re- 
spect to the hopping integral. In the zero order (t = 0) 
the electron is localized either on the left or on the right 
site, so that 



a'i(« m ) = exp 



E 



/m(l) 

Muj 2 



,a 2 (u m ) = 
(7) 



or 



a-Kum) = 0, a r 2 (u m ) = cxp 



Mto^f f m {2) 



(8) 

In the first order a±(u m ) and a 2 (u m ) are linear combi- 
nations of aj(ii m ) and a5(u m ) as 



ai{u m ) 



ol( o"™ ) ) + "(«: 
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(9) 
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FIG. 2: Schematic representation of the Frohlich (or Extended 
Holstein) model and the Holstein model. The electron on the 
site 1 interacts (a) with the sites m = — 1 and m = and (b) 
with only site m = 1 of the ion chain, in the Frohlich and the 
Holstein model, respecively. 



Substituting Eq.(9) into Eq.(5) and Eq.(6) one gets a 
system of linear equations with respect to a and j3. The 
standard secular equation for the energy is 



det 



E — Nuj/2 + E p 



t 



t 



E - Noj/2 + E p 



= (10) 



where E p = J2 m f^ l (l)/2Mu; 2 is the polaronic shift, N 
is the number of ions in the chain and 



t = t 



o' 1 (u m )a2(w m )dw m / / \a[(u m )\ 2 du 



(11) 



is a renormalised hopping integral. Then the lowest en- 
ergy levels of the system are found as E± — Nui/2—E p ±t. 
The evaluation of Eq.(ll) with explicit a l x (u m ) and 
a 2 {u m ) results in t = £exp (— g 2 ), where 
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2Muj 3 ^ 



- / m (l)/ m (2)) 



(12) 



The same expression for polaronic shift was obtained 
in Ref. |33j| by using the canonical Lang-Firsov transfor- 
mation. We would like to stress that the model yields 
less renormalization of the effective mass than the Hol- 
stein model. For simplicity let us take into account only 
nearest-neighbors interactions, as Fig. 2a. In this case our 
model yields E p = /g (1)/Muj 2 and the mass renormaliza- 
tion m* I'm = exp(E p /2oj), while the Holstein model with 
the local interaction, Fig. 26, yields m* jm = exp(E p /oj). 

The factor 1/2 in the exponent provides much lighter 
Frohlich small polarons compared with the Holstein 
model, Fig. 3. If one considers the Coulomb-like in- 
teraction with the whole upper chain, one gets the factor 
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FIG. 3: Ratio of masses of nonadiabatic Frohlich and Holstein 
small polarons as a function of the electron-phonon coupling 
constant A at different values of uj/t. 



0.28 j^H instead of 0.5 in the exponent, that means even 
less renormalised effective mass. 

In the opposite adiabatic regime we use the Born- 
Oppenheimer approximation representing the wave func- 
tion as a product of wave functions describing the "vi- 
brating" ions, x( u m)i an d the electron with a "frozen" 

ion displacements, ( i/j(u m ) f{u m ) ) (T means trans- 
pose matrix). Terms with the first and second derivatives 
of the "electronic" functions ip(u m ) and ip(u m ) are small 
compared with the corresponding terms with derivatives 
of x(u m )- The wave function of the "frozen" state obeys 
the following equations 



E(u m ) 



tip{u m ) + 



', /m(l)ltr, 



ip(u m ) - tcp(u m ) = (13) 



E(u m )-J2f m (2) 



<p(u m ) = 0. (14) 



The lowest energy is 

E{u m ) = 1^/ + Um - [J(£/m«m) 2 +* 2 ] 1/2 , (15) 

m m 

that plays a role of potential energy in the equation for 

X(«m) 



X (u m ) = 0. (16) 



Here /+ = [/ m (l) + / m (2)] and /" = 
[/m(l) — /m(2)]. Because of the infinite number of 
variables, Eq.(16) can not be reduced to a double-well 
potential problem However, in the nearest-neighbor 
approximation with three relevant ions in the upper 
row, Fig. 2a, it can. Indeed, using the transformation 
£,m = u m + u -m and rj m = u m — U— m and = ,f^ m , 
f~ = — /~ j one can rewrite the latter equation as 



3 d 2 
E-2lu+ -E p + 



u(vi) x(m) = o. (17) 



Here 



U(vi) 



2 9 

^ m 



1/2 



(18) 



is the familiar double-well potential, and /i = M/2. Then 
the determination of the energy splitting is similar to 
the case considered in The splitting appears due to 
penetrable character of the double-well potential: 



AE 



■ exp 



'ftp 

-2 J \ P {m)\dm 

o 



(19) 



where r\t v = 1 — 1/ \J fiu> is the classical turning point 
corresponding to the energy of the ground state Eq = 
uj/2 + U min (r)i) and p(r}i) = y/2^i(E - Ufa)) is the 
classical momentum. The result is AE — Aexp(— g 2 F ) : 
where 



7T 



2u 



1 - 



2u 



,3/2 



(20) 



and 



9f 



E p 
2uj' 



2 = ±P K l/2 1 



2u 



(21) 



Here uj = lu^/k is the renormalised phonon frequency, 
k = (1 - 1/A 2 ) and A = E p /(2t). In the Holstein model 
one should replace E p in Eq.(20) and Eq.(21) by 2E p . 
The comparison of two models shows that the polaron 
in the Frohlich model remains much lighter, than the 
Holstein polaron also in adiabatic regime, Fig. 4. 

In conclusion, we have solved an extended Holstein 
model with a long-range Frohlich interaction in the adia- 
batic limit. We have found the hopping integral of SFP, 
and compared it with SHP. The small adiabatic Frohlich 
polaron is found many orders of magnitude lighter than 
the small Holstein polaron both in the nonadiabatic and 
adiabatic, Fig. 4, regimes. One of us (B.Ya.Ya) is grateful 
to NATO and the Royal Society for their financial sup- 
port (grant PHJ-T3). We also acknowledge the support 
of Lever hulme Trust, UK. 
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FIG. 4: Ratio of the SFP mass to the SHP mass as a function 
of A in the adiabatic regime. 
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